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Abstract 

In this work we derive the copulas related to Manneville-Pomeau processes. We examine 
both bidimensional and multidimensional cases and derive some properties for the related 
^-H copulas. Computational issues, approximations and random variate generation problems are 

also addressed and simple numerical experiments to test the approximations developed are 
also perform. In particular, we propose an approximation to the copula which we show to 
^ converge uniformly to the true copula. To illustrate the usefulness of the theory, we derive 

^ a fast procedure to estimate the underlying parameter in Manneville-Pomeau processes. 

QQ Keywords. Copulas; Manneville-Pomeau Processes; Invariant Measures; Parametric Esti- 

y—( mation. 

H 

C/^ 1 Introduction 

Cd The statistics of stochastic processes derived from dynamical systems has seen a grown attention 

^ in the last decade or so (see Chazottes et al. (2005) and references therein). The relationship 

between copulas and areas such ergodic theory and dynamical systems also have seen some 
fSl development, especially in the last few years (see, for instance, Kolesarova et al. (2008)). In 

this work our aim is to contribute with the area by identifying and studying the copulas related 
to random vectors coming from the so-called Manneville-Pomeau processes, which are obtained 
as iterations of the Manneville-Pomeau transformation to a specific chosen random variable (see 
Definitions 2.1 and 2.2). We cover both, bidimensional and n-dimensional cases, which share a 



lot more in common than one could expect. 

The copulas derived here depend on a probability measure which has no closed formula. In 
order to minimize this deficiency, we propose an approximation to the copula which we show 
to converge uniformly to the true copula. The copula also depend on several functions which 
^ have to be approximated as well, so the approximation depends on several intermediate steps. 

^ The results related to the convergence of the proposed approximation presented here are far 

more general than we need and actually allows one to change these intermediate approximations 
and still obtain the uniform convergence result for the approximated copula. We also address 
problems related to random variate generation of the copula and present the results of some 
simple numerical experiments in order to assess the stability and precision of the intermediate 
approximations. The usefulness of the theory is illustrated by a simple application to the problem 
of estimating the underlying parameter in Manneville-Pomeau processes. 

The paper is organized as follows: in the next section, we briefly review some concepts 
and results on Manneville-Pomeau transformations and processes and on copulas. Section 3 
is devoted to determine the copulas related to any pair {Xt,Xt+h) from a Manneville-Pomeau 
process and to explore some consequences. In Section 4, the multidimensional extensions are 
shown. In Section 5 an approximation to the copulas derived in Section 3 is proposed. This 
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approximation, which is shown to converge uniformly to the true copula, is then applied to exploit 
some characteristics of the copulas related to Manneville-Pomeau process through statistical 
and graphical analysis. Some computational and random variate generation problems are also 
addressed. In Section 6 we illustrate the usefulness of the theory by deriving a fast procedure to 
estimate the underlying parameter in Manneville-Pomeau processes. Conclusions are reserved 
to Section 7. 



2 Some Background 

In this section we shall briefly review some basic results on Manneville-Pomeau transformations 
and related processes as well as some concepts on copulas needed later. We start with the 
definition of the Manneville-Pomeau transformation. 

Definition 2.1. The map : [0, 1] — > [0, 1], given by 

T,{x) = aj + x^+'Cmodl), 
for s > 0, is called the Manneville-Pomeau transformation (MP transformation, for short). 

In what follows, A shall denote the Lebesgue measure in I := [0, 1] and the A;-fold compo- 
sition will be denoted, as usual, hy = o • • • o T^. Figure 1 shows the plot of the MP 
transformation for the values of s G {0.5, 1, 10, 100}. The plots show the usual behavior of the 
MP transformations: for any s, they are increasing and differentiable functions by parts in /. 
Furthermore, for any s > 0, the function will have exactly 2^^ parts. 







Figure 2.1: Plot of the Manneville-Pomeau transformation for different values of s G {0.5, 1, 10, 100}. 



Pianigiani (1980) shows the existence of a T^-invariant and absolutely continuous measure 
with respect to the Lebesgue measure in / which will be denoted henceforth by /u^. However, 
the proof uses Perron-Frobenius operator theory and is, for practical purposes, non-constructive 
so that an explicit form for a T^-invariant measure is unknown. However, this measure will be 
a Sinai-Bowen-Ruelle (SBR) measure in the sense that the weak convergence 



^ n—l 



(2.1) 



k=0 



holds for almost all x £ / and all ^g-continuity set^^, where 6a{-) is the Dirac measure at a. 

As a dynamical system, the triple (/, /is,Ts) is exact (that is, limfc_j.oo (/^s ° Tg){A) = 1, for 
all positive //^-measurable sets A) which implies ergodicity and strong- mixing. When s < 1, 



^Recall that a set A is a /^-continuity set if fi{dA) — 0, where dA denotes the boundary of A. The measure 
theoretical results applied here can be found, for instance, in Royden (1988). A good reference in weak convergence 
of probability measures is Billingsley (1999) and for ergodic theoretical related results, see PoUicott and Yuri 
(1998). 
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lis is a probability measure, while if s > 1, //^ is no longer finite, but cr-finite (see Fisher 
and Lopes (2001)). Furthermore, it can be shown that /i^ has a positive, bounded continuous 
Radon-Nikodym derivative d^^ = hs(x)dx, fact that will be useful later. For further details 
in the theory of MP transformations and related results, we refer to Pianigiani (1980), Young 
(1999), Maes et al. (2000) and Fisher and Lopes (2001). For appUcations, see Zebrowsky (2001), 
Olbermann et al. (2007) and Lopes and Lopes (1998). 

Definition 2.2. Let s G (0, 1) and let Uq be a random variable distributed according to (the 
probability measure) fig- Let (p : [0, 1] — > E, be a function in £.^{ns). The stochastic process 
given by 

Xt = {ipoT*){UQ), for all i e IN, 
is called a Manneville-Pomeau process (or MP process, for short) . 

The MP process, as defined above, is stationary since fig is a T^-invariant measure and 
/Is ^ A. It is also ergodic since /is is ergodic for Tg. By its turn, copulas are distribution 
functions whose marginals are uniformly distributed on /. The copula literature has grown 
enormously in the last decade, especially in terms of empirical applications and have become 
standard tools in financial data analysis (see Nelsen (2006) and references therein). The next 
theorem, known as Sklar's theorem, is the key result for copulas and elucidates the role played 
by them. See Schweizer and Sklar (2005) for a proof. 

Theorem 2.1 (Sklar). Let Xi, ■ ■ ■ ,X„ be random variables with marginals Fi, ■ ■ ■ ,Fn, respec- 
tively, and joint distribution function H. Then, there exists a copula C such that, 

H{xi,- ■ ■ ,Xn) = C{Fi{xi),- ■ ■ ,Fn{xn)), for all (xi, • • • ,x„) G R"'. 

If the Fi 's are continuous, then C is unique. Otherwise, C is uniquely determined on Ran(Fi) x 
• • • X Ran(F„). The converse also holds. Furthermore, 

C{uu--- ,Un)=H{F!;-^\ui),--- ,F^-^\un)), for all (m,--- ,n„) G/", 

where for a function F , F^~^^ denotes its pseudo-inverse given by F^~^\x) := inf \u G Ran(F) : 
F{u) > x]. 

The next theorem, whose proof can be found, for instance, throughout Nelsen (2006), shall 
prove very useful in what follows. Except stated otherwise, the measure implicit to phrases like 
"almost sure", "almost everywhere" and so on will be the (appropriate) Lebcsgue measure. 

Theorem 2.2. Let X and Y be continuous random variables with copula C. If f is an almost 
everywhere decreasing function then C f(^x),Y{'<^-,v) = u — Cx,y{u, 1 — v). Furthermore, if fi and 
/2 are functions increasing almost everywhere, then C f^(^x)j2{Y){u-,v) = Cx,y{u,v) . 

For an introduction to copulas, we refer the reader to Nelsen (2006). For more details and 
extensions to the multivariate case with emphasis in modeling and dependence concepts, see 
Joe (1997). The theory of copulas is also intimately related to the theory of probabilistic metric 
spaces, see Schweizer and Sklar (2005) for more details in this matter. 

3 Copulas and MP Processes: Bidimensional Case 

In this section we shall investigate the bidimensional copulas associated to pairs of random 
variables coming from MP processes which we shall call MP copulas. As we will see later, the 
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multidimensional case is very similar to the bidimensional case, so we shall give special attention 
to the latter. 

First, let {Xn}neTN be an MP process with parameter s G (0,1) and cp G ^^{fJ-s) be an 
increasing almost everywhere function. Throughout this section and in the rest of the paper, we 
shall treat s G (0, 1) as a given fixed number. Let 

Fo{x) :=P(C/o <a:) =Ms([0,a;]). 

Since fis ^ A, /x^ is non-atomic and, therefore, Fq will be (uniformly) continuous. The existence 
of a positive Radon-Nikodym density for also shows that Fq will be increasing and its inverse 
will be well defined. Let Ft be the distribution function of T*(C/o), for all t G M. For x £ I, 
notice that 

Ftix) :-P(Ti((7o)<x) =/i,((T*)-i([0,x-])) =M.([0,x]) =Fo(x), (3.1) 

since is a Ts-invariant measure. 

In what follows, we shall need the solution for the inequality Tg{X) < y, y £ (0, 1), in X, 
for X a random variable taking values in /. Now, since each of the 2* parts of T* is one-to-one 
in its domain, the inverse of T* will also be continuous by parts and each part will also be a 
one-to-one function in its domain. Let = atfi, • • • , 2* = 1 be the end points of each part of 
T*. We shall call each interval [at^fc) a node of Tj, for k = 0, - ■ ■ , 2* — 1 and t > 0. The 

(piecewise) inverse of T* can be conveniently written as 

y ^ {TtAy)r-- ,Tt.2^~i{y)), (3.2) 

where Tt^y) denotes the inverse of T* restricted to its fc-th node, for all A; G {0, • • • , 2* — 1}. 
Notice that both Tt j^ and at^k depend on s for each /c, but since no confusion will arise, and for 
the sake of simplicity, we shall omit this dependence from the notation as we shall do in several 
other occasions. Now, the solution of the inequality Tl{X) < y \n X can be determined and is 
given by X G Atfl{y) U • • • U A,2^-i{y), where 

At,k{y) =[at,k,Tt,k{y)\, (3.3) 

which will be a proper closed subinterval of [at,fc, at,fc+i), for each A; = 0, • • • , 2* — 1. Notice that 
At,k{y) (whose dependence on s was omitted from the notation) is just the inverse image of 
[0, y] by the transformation T* restricted to the node [at^ki0.t,k+i)- We can now use this result 
to prove the following useful lemma. 

Lemma 3.1. Let X he a random variable taking values in I and let Tg he the MP transformation 
with parameter s > 0. Then, for any t G IN and x G I, 

2'-l 
k=0 



where At^k's are given hy (3.3). 



Proof: The result follows easily from what was just discussed and from the fact that the 
intervals At^kS are (pairwise) disjoints. ■ 

As for the copulas related to MP processes, in view of the stationarity of the MP process, 
the following result follows easily. 

Proposition 3.1. Let {Xn}neJM he an MP process with parameter s G (0, 1) and ip G C^{fis) be 
an almost everywhere increasing function. Then, for any t, /i G M, 

Cxt,Xt+Hiu,v) = Cxo,xdu,v), 

everywhere in . 
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Proof: As consequence of the stationarity of {XjjtgiN, if we let the joint distribution of the 
pair {Xp, Xg) for any p, q G W, p / g, be denoted by Hp^g{-, •), it follows that for all x,y £ (0, 1), 
t G I N an d ^ G IN* := ]N\{0}, Ht^t+h{x,y) = HQ^h{x-,y). Now, upon applying Sklar's Theorem 



and (3.1), it follows that 

for all {u,v) e P. ■ 

Corollary 3.1. Let Tg he the MP transformation for some s G (0, 1), /Ug he a Tg-invariant 
prohahility measure and let Uq he distributed as fis- Then, for any t, /i G IN, /i / 0, 



a 

everywhere in . 



T*((7o),r*+''((7o)(^'^) - C'^o,Th(t/o)(^'^) 



Proof: Immediate from Theorem 2.2 applied to Proposition 3.1 



Now we turn our attention to determine the copula associated to any pair of random variables 
{Xp,Xq), p,q G¥i, obtained from an MP process with (p increasing almost everywhere. For the 
sake of simplicity, let us introduce the following functions: let h he a positive integer and for 
k = (),■■■ , 2^* - 1, let Th,k ■■ I [Fo{ah,k), -Fo(aft./c+i)]be given by 

Th,k{^) :=Fo{Th.k{F^\x))). 

Notice that for each k, J^h,kiO) = -fb(fl/i,fe) and J^h,ki^) = -^o(«/i,A:+i) and J^h^k is a one to one, 
increasing and uniformly continuous function. 

Proposition 3.2. Let {Xn}neJM be an MP process with parameter s G (0, 1), G C^{^s) he an 
increasing almost everywhere function and let Fq be the distribution function oJUq. Then, for 
any t, /i G IN, h ^ and (u, v) G I^, 

/no-l \ 

Cxt,Xt+h(u,w) = I ^ Fh,k{v) - i^o(a/i.fe)|'^iN-('^o) + min{w, J'/j^nJt;)} -Fo(a/j,„J, (3.4) 

\ k=0 ) 

where equals 1, if x ^ IN* and 0, otherwise, {ah,k}'kLQ are the end points of the nodes of 

T^ and no := no{u;h) = {k : u e [Fo{ah^k), Fo{ah,k+i))} G {0, • • • , 2^* - 1}. 



Proof: By Propositions 3.1 and 2.2 it suffices to derive the copula of the pair (^Uo,T^{Uo)) ■ 
So let again {Xn}n&f^ be an MP process with parameter s G (0,1) and if G C^{fJ.s) be an 
increasing almost everywhere function and let ifo,/t(', •) denote the distribution function of the 
pair {Uo,T^{Uo)). Notice that 

HoAx,y) =P(C/o < x,T^{Uo) < y) = P(C/o < x,Uo G UllV ^hMv)) 

= P(C/oe [O,x]nuto'^'«,fe(y)) =P(t^oeUto'[[0'^]nA,fe(l/)]) 

= p(c/oG [o,.T]nA,fc(2/)), 

fc=0 

for any x,y G (0,1). Now let ni := ni{x;h) = {A; : x G [ah,k,Qh,k+i)} £ {0, • • • ,2^ — 1} and 
assume for the moment that ni > 1. Since Ah^kiu) = [o'h,k,Th,k{y)], it follows 

ni — 1 
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ni — 1 



= ^ Hs{Ah.k{y)) + f^s{[ah,ni,Th,nAy)]r\[ah,ni,x]) 
fe=0 
ni — 1 

= X! fJ's{[ah^k,Th.k{y)]) + ^J.s{[ah.nl,n-lm{x,Th^nl{y)}])l 

fc=0 

which can be written, since -^0(2^) = /f^s([0,2;]) is increasing, as 

ni — 1 

Ho,h{x,y) = ^ [Fo{Th,k{y)) - Fo{ah,k)] + min {Fo(x), Fo(7A,„j (y)) } -Fo(a;,_„J. 
fe=o 

If ni = 0, the summation is absent of the formula and we have 

Ha,h{x,y) = v[)\n{Fo{x),Fo(Th.,o{y))] ~ ^0(0^,0), 
so that, in any case, we have 

Ha,h{x,y)^ ( ^ [i^o(7^,fc(y)) - Fo(a,i_iO] ]^iN*("i) +min{Fo(a;),i^o(7^,ni(y))} - ^o(a/i,ni)- 
V fc=o / 

Now upon applying Sklar's Theorem, it follows that 

Cu.,niUo){n,v) = H^jXFo\^IFh\^)) = H^,h{F^\u),F^\v)) 

/no-l \ 



fc=0 



where no := no{u;h) = ni(FQ = {A; : n G [FQ{ah,k), Fo{o.h,k+i))} ■ The result now 

follows from Proposition |3.1[ ■ 



Remark 3.1. Notice that the copula (3.4) can be expressed in terms of as 



C 



'xuXt+du,v) = i^^^s(^[ahM,Th.k{FfJ^{v))]'j\^JMinQ) + 



. k=0 



+ A's([aft,no,niiii{i^o ^{u),Th,no{F„ (3.5) 



which will prove useful in Section 5. Also, expression (3.5) is helpful if one desires to verify 



directly that the marginals of (3.4) are indeed uniform. 



In the next proposition we address the case where ip is an almost everywhere decreasing 
function. In view of Theorem 2.2, one could, at first glance, think that a result like Cxo,x,i = 



Cxt,Xt+h would not hold anymore, but in fact it still does, as it is shown in the next proposition. 

Proposition 3.3. Let {XnjneiN an MP process with parameter s G (0,1), if £ ^^(fJ-s) be 
an almost everywhere decreasing function and let Fq be the distribution function ofUo. Then, 
Cxo,Xh('U5 "w) = C'x(,Xt+h(^^; ^) everywhere in P and, for any /i G IN and h ^ 0, 

Cxt,Xt+hi.'^^^) = U + V-1+ [-^fe,fc(l -v)- Fa{ah,k)] ^w(?to) + 

\A:=0 / 

+ minjl - u,Fh,na{^ - ~ Fo{ah,no), (3.6) 

for all {u,v) G , where {cLh,k}'k=Q ^he end points of the nodes of and uq := no{u; h) = 
{k:ue{l- Fo(a;i_fc+i), 1 - -^o(a/i,fe)] } • 
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Proof: Since the inverse of an almost everywhere decrea sing function is still decreasing almost 



everywhere and Xt = ip(Tg{Uo)) , upon applying Theorem 2.2 twice, it follows that 

= u- {l-v-Cxt,Xt+,A^-u,l-v)), 
or, equivalently (changing n by 1 — -u and f by 1 — v), 

CxuXt+du, v)^u + v-l + Cj,t^jj^^,j,t+H^jj^^{l -u,l-v). (3.7) 

iting equation 
and Theorem 



Now (3.6) follo ws upon applying Proposition 3.2 with the identity map and sub stitu ting equation 
(3.4) into (3.7). As for the equality Cx^.x^ v^i ^) = C'xt,Xt+h('U5 v). Corollary 



2.2 applied to (3.7) yield 



3.1 



Cxt^Xt+Hiu,-") = u + i;-l + Cc/„^Ti'((7o)(l-w>l-w) 

= U + V -1 + C^-i(^(t/o)),v-i(¥'(T^''((7o)))(l - U,l-v) 

everywhere in I^, as desired. ■ 

Remark 3.2. In view of the "stationarity" results of Theorems |3.1| and |3.3[ a copula associated 
to a pair {Xt^Xfj^h) from an MP process will be referred as lag h MP copula. 



The copulas in (3.4) and (3.6) are both singular, as it can be readily verified. So the question 
that naturally arises is, for each h, what is the support of Cxt,Xt+h^ The question is addressed 
in the next proposition, which will be useful in Sections 5 and 6. For simplicity, for a given MP 
process and > 0, let ^^^j^^^ : [-^o(o/i,/c), -Fo(o/i,/c+i)) — >• / be functions defined by 

for all /c = 0, • • • ,2^ — 1. Notice that, for each k, is the linear function connecting the points 
{Foiah,k),0) and (Fo(oh,fc+i), l), while connects the points {Fo{ah,k), l) and {Fo{ah,k+i), O). 

Proposition 3.4. Let {Xn}ne'M be an MP process with parameter s £ (0,1), for ^pi £ C^ifis) 
an almost everywhere increasing function and let {l^jnew be an MP process with parameter s G 
(0, 1), for S >C^(/is) o-n almost everywhere decreasing function. Also let Fq be the distribution 
function of Uq. Then, for any t,h h > 0, 

snpp{Cx„x,+J = U to'{("'<fe(^^)) ■■ " e [Foian,k),Fo{an,k+i))} (3.8) 

and 

supv{Cy,,y,^J = Uf=oHK^^,feH) ■■ " G [Foiah.k),Fo{ah,k+i))}- (3.9) 

Proof: Let R = [ui,U2] x [^1,^2] be a rectangle in P and let its Cxt,Xt+h-volume be denoted 
by {R)- Let A; G {0, • • • , 2^* — 1} be fixed and suppose that Ui G [-^o(o/i,fc), -Fb(«?i,fc+i)] ■ This 
implies that uq = k for all four terms in Vcx-{R), hence the summands and constants on the 
copula cancel out so that we have 

y^CxiR) = inin{wi, J'^^fe(ui)} + min {u2, J^h,k{v2)} - min {ui, J"ft_fe(w2)} - min {7/2, -7^h,A;(wi)} 

where M{u, v) = min{u, v} is the Frechet upper bound copula whose support is the main 
diagonal in P. Since [ui,U2] x [Fh,k{vi), Fh,k{vi)] C [i^o(afe,fc), ^o(a^fc+i)]^, ^CxiF.) > if, and 
only if, Rf]{{u,ilf^{u)) : u e [Fo{ah,k), Fo{ah,k+i))} ^ ^- 



8 



Copulas Related to Manneville-Pomeau Processes 



Analogously, denoting the Cyj^yj^^-volume of R by Vcy{R), if G [l — Fo{ah,k),^ — 
Fo{ah,k+i)], we have 

Vcy{R) = min{l - ui,J'h,k{^ - ^'i)} + min {l - U2,-F/j,/c(l 

- min {1 - ui,Th^k{l - W2)} - niinjl - W2,-^/i,fc(l - Vi)} 
= Vm{[1 - Ml, 1 - U2] X - W2), J"h.fe(l - «i)]). (3.10) 

Since [1 - ni,l - U2] x [J)i,a..(1 - - V2)] C [Fo(afe,fc), i^o(a/i,fc+i)]^, is pos- 

itive if, and only if, R£\{[u,iJ^i^{u)) : u G [-Fo(a/i,fc), -Po(ah,fc+i)) } / (notice the terms 
1 — Vi in expression ( |3.10[ ), for i = 1,2). Now (3.8) and ( |3.9| ) follow by observing that / = 

ULo^ [Fo{ah,k),Fo{ah,k+i)] = ULo^ ~ Fo{0'h,k+i), 1 - ^o(a?i,fc)] • ■ 
Remark 3.3. We end up this section by noticing that as an application of Propositions |3.1| 



and 3.3, together with the so-called copula version of Hoeffding's lemma (see Nelsen (2006)), 
we can show in a rather different way that an MP process is weakly stationary. Let Fx^ be the 
distribution function of Xt and notice that Fxt{x) = Fxq{x), for all t G W, by the stationarity 
of {Xt}tejM and since Cxt,Xt+h ("U, w) = Cxo,Xh(^, "w), the result follows immediately. 



4 Multidimensional Case 

In this section we are interested in extending the results from the previous section to the mul- 
tidimensional case, that is, in this section we are interested in deriving the copulas associated 
to n-dimensional vectors (Xt^,--- ,Xt„), ti,--- ,tn £ IN, coming from an MP process with if 



an increasing almost everywhere function. In view of Theorem 2.2, it suffices to derive the 
copula associated to the vector (Tg^{Uo), ■ ■ ■ ,T*"(J7o))- It turns out that there are more simi- 
larities between the bidimensional and multidimensional cases than one could expect. In fact. 



an expression very similar in form to (3.4) holds for the multidimensional case as well. 

Let {Xn}n£K be an MP process with parameter s G (0,1) and f S C^ins) be an al- 
most everywhere increasing function. For the sake of simplicity, we shall use the following 
notation: let a, 6 G M, a < 6, we shall write Xa-h '■= (x^,--- ,Xf,) and for a function /, 
f{xa:b) '■= {f{xa), • ' ' , fixb))- Again we shall denote the distribution function of Uq by Fq. 

Theorem 4.1. Let {Xn}n£K be an MP process with parameter s G (0,1), with (/? G C^{fis) o-n 
almost everywhere increasing function. Let ti, • • • ,in G IN and set hi := ti — ti. Then, for all 
{ui,--- ,Un) G 



/' no~l 



U^,- ^Xt^ui, - ■ ■ ,u„)= ( ^ FQ{hh,„k{FQ \u2:«))) - -Fo(aft„,fe)UiN'("-o) + 



fc=o 



ah„,nQ ), (4.1) 



where uq := no{ui,n) = {k : m e [Fo(a/i„,fc), -^o(ah„,fe+i)) }, {ah„^k}f=o are the end points of 
the nodes of T^"^, for i = 2, • • • ,n, j = 0, • • • 2^' — 1, Fh^j is given by (3.2) and for a vector 



(x2, • • • ,x„) G /" ^ , bh^,k{x2:n) = min {ci{xi;hn,k)], with 



{^i , hji , /c) 



a/i„,fe, if Bi{xi;hn,k) 

Bi{xi;h„,k), otherwise. 



and 



Bi{xt;h„,k) = min {Fhi,j{xi) ■ Fhijixi) > ah„,k and ah,,j < ah„,k+i}- 

j=0,--- ,2'^i-l 
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Proof: Let {Xn}nem be an MP process with parameter s G (0, 1) and (p G C^^fji-s) be an almost 
everywhere increasi ng fu nction. Without loss of generality, we can assume that < ti < • • • < t^- 
In view of Theorem 2.2, it suffices to work with the vector (TI^{Uq), ■ ■ ■ ,r*"([/o)). Let Htj^^...^t„ 
be the distribution function of (Tg^{Uo), ■ • • ,T*"(f7o))- Let hi = U — ti, for each i = 1, ■ ■ ■ ,n, 
and notice that hi > since ti < ti, for alH = 2, • • • , n. Let (xi, • • • , x„) G (0, 1)" and for the 
sake of simplicity, let Yt-^ := T^^{Uo), so that we have 

= p(rt,e [0,a;i],yt,e ur=V'^/..,fc(a;2), • • • ur=T'^h„,fc(x„)) 
= p(rt,e [0, xiin r=2 [U tlo^An^.kix,)] ) 

= p(c/o e nr=2UfelV'[[0'^i]n^/>„fc(a:0]), (4.2) 



where A^^^k^ are given by (3.3) and the last equality is a consequence of the Ts-invariance of 
Us- For A; = 0, ••• ,2^"-^, let 

Ah„,k{.X2:n) = A„.fc(a;„)nr=2MU|=o'^^/».j(a;i)]- 
In order to simplify the notation, for i = 2, • • • , n and A; = 0, • • • , 2'^" — 1, let 

Bi{xi; hn, k) = min {Th,.j[xi) : Th-.jixi) > at^^k and ah„j < ah„.k+i}- 

j=0,-- ,2''i-l ■ 

For each k and i, Bi(xi;hn,k) is either the smallest Thi,j{xi) which is greater than ah„,k and 
such that the correspondent Af^j^Xi) has non-empty intersection with Ah„^k{xn), or empty. Let 



Ci (^Xi , /i^i , /l) — 



a/i„,fe, if B,{xi;hn,k) 

Bi{xi;h„,k), otherwise. 



Then, for each k = 1, - ■ ■ , 2^" — 1, setting bh„ k{x2:n) = min {ci{xi; hn, k)\, it follows that 

' i=2,--- ,n 

Ah„,k{x2:n) = [ah„,k,bh„.k{x2:n)], 

which is a closed subset of [ah„,k, a/i„,fc+i]- Also notice that, from the definition of bh„,kix2:n), we 

could have Ah„^kix2:n) = {a/i„,fc}, in which case we set ^h„,fc(aJ2:n) = (although from a measure- 
theoretical point of view, this correction makes no difference). Again we are omitting the 

dependence in s from the notation on both, bh^^k and Ah„^k- Each bh„,k{x2:n) above determines 
the smallest Thi,j{xi) that lies on the fe-th node of T^^^ (which has the smallest nodes among all 
Tg^^s), so that Ah^^s are just the intersection of all Ah^^k{xi)'s with end point in the k-th node 
of Tg". Also notice that the A^^^kS are pairwise disjoints. One can rewrite (4.2) as 

Ft,,.. ,t„ (a^i , • • • , x„) = P (c/o e U 'IV [Ah,Mx2:n)r\ [0, ^i]] ) . (4.3) 

Now, let ni := ni{xi;n) = \k : xi G [ah„,k,o-hn,k+i)} G {0, • • • ,2^" — 1}, and assume for the 
moment that rii > 1. Then (4.3) becomes 

ni — 1 

-ffii,-,t„(a;i, ■ ■ • = ^ P(C/o e ^/i„,fe(2;2:n)) +P(C/o e A,i^^„,(x2:„)n [ah„,«i,a;i]) 

fc=0 
ni — 1 

= ^ Ms(K„,fc,&/i„,fc(a;2:n)]) +Ms([a/i„,™i,min{a;i,6/,„^„j(a;2:„)}]) 

/c=0 
ni — 1 

= X] [^o{bh„,k{x2:n)) - Fo{ah„,k)] +min{Fo(a;i),i^o(^;i„,«i(a;2:n))} - ^o(a/i„,ni)- 



/c=0 
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If rii = 0, then 

Ht^,... .t„{xi, ■ ■ ■ ,Xn) = mm {Fo{xi), FQ{bh^,Q{x2:n))} - ^o(a/i,.,o)- 
In any case, we can write 



. fc=0 



+ mm {Fo{xi), Fa{bh„^n^{x2:n))} - -Fo(aft„,„i). 

Recall that the distribution function of T*(C/o) is also Fq by the Tg-invariance of /i^. Now 
applying Sklar's Theorem, it follows that, 



Cxt,,-,Xt^iui,--- ,Un) ^ Ht,^...,t„{F^^{ui),--- ,F^^{Un)) 

+ min{ui,Fo(6/i„,„(,(i^(7^(u2:n)))} - Fo(aft„,„o)- 
where no := ni{F^^^{ui),n) = {k : m £ [Fo{ah„,k), Foio-hr^k+i))} , which is the desired formula. 



Remark 4.1. Notice that the proof of Theorem 4.1 from equation (4.3) on is exactly the same 



as the one in Proposition |3.2| with the obvious notational adaptations. 

Now we turn our attention to the case where ip is an almost everywhere decreasing function. 



In view of Theorem 2.2, one cannot expect a simple expression for the copula. What happens 
is that the copula in this case will be the sum of the lower dimensions copulas related to the 
iterations Tg{Uo), as the next proposition shows. 

Proposition 4.1. Let {Xn}n£TN be an MP process with parameter s G (0, 1), and tp £ C^{fis) 
he an almost everywhere decreasing function. Let t, hi, ■ ■ ■ , hn G M, < hi < ••• < /i„, 
and set Yq := Uq and := T^*([/o)- Denote the copula associated to the random vector 
{Xt, Xt+hi, • • • jXt+hn) by Ct- Then the following relation holds 



n n 



Ct{uo, • • • , Un) = l-n + ^Ui + ^ ^ C:Yi,Yj (1 - Ui, 1 

i=0 1=0 j=i+l 



n n 



fci=0fc2=fcl+l /£„_l=fc„_2 + l 

+ (-irCy„,y,,..,y„(l-UO,--- (4.4) 



everywhere in 

Proof: Let t, /ii, • • • , /i„ G M, < /ii < • • • < /i„, t / 0. Set Yq := Uq, Ifc := T^^iUo) and 
yk ■■= ^{xk)- We have 

) = P(C/o > y<d,Yl > Vl, - ■ ■ ,Yn > Vn) 

= P(c/o > yo\Yi > yi,Y2 >y2,--- ,Yn> y„)p(Yi > yi, • • • , r„ > zj„) 

= P(Yi>yi,--- ,Yn>yn) -I'{Uo<yo,Yi>yi,--- ,Yn>yn). (4.5) 



Upon applying a long chain of a conditioning argument on both terms in (4.5), we arrive at 



Hxo,X^^,-,X^^{xo,Xi,--- ,Xn) = l-^Fo{yi)+Y^ ^ HY^,Y,{y^■,y]) + 

1=0 i=0 j=i+l 
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fel=0fc2 = fel+l /c„_i=fe„_2 + l 

+ (-l)"i/,7o,yi,-,y„(l/o,--- 



(4.6) 



A simple calculation (using the T^-invariance of fig) shows that, for all tSlN* and x G (0, 1) 

FxAx) = I - Fa{^{x)) and F^^ (x) ^ {F^\l - x)) , 
so that, the result follows upon applying Sklar's Theorem to (|4.6|) (recall that yk = ip{xk)). 



Remark 4.2. Notice that the copula in Proposition 4.1 can be explicitly calculated since (4.4) 
is written as sums of the copulas of vectors containing Uq and T*([/o) for different t's, so that 



the desired formulas can be deduced in terms of the copulas in Theorem 4.1 



5 Numerical Approximations to the MP copulas 

The MP copulas derived in the last sections do not have readily computable formulas, especially 
because fis does not have explicit expression and because even apparently simple tasks like 
determining the discontinuity points of or to compute explicit formulas for the branches 
of Tg can be highly complex ones. However, one can still study the copulas derived in the 
last sections by using appropriate approximations to the functions appearing in the copula 
expression. Besides the invariant measure /x^, computation of the bidimensional copulas so far 
discussed also involves computation of the quantile function Fq^, the inverse of and the end 
points {a/i,fc}fc=o nodes of T^. 

In this section our goal is to derive simple approximations to these functions in order to 
obtain an approximation to the copula itself, which we shall prove to converge uniformly in its 
arguments to the true copula. The approximations presented here are simple ones, usually a 
linear interpolation based on a grid of values, but the technique and results we shall use and 
prove here are stronger and cover a wide range of approximations, for instance, all results hold if 
we use some type of spline interpolation instead of a linear one. This is so because the functions 
to be approximated are generally very smooth. We also evaluate the stability and performance 
of the approximations by simple numerical experiments. 



Approximation to /i^ 

We start with an approximation to fig- In this direction there are at least two ways to compute 
approximations to fig- One way is by using the ideas and results outlined in Dellnitz and Junge 
(1999), which are based on a discretization of the Perron- Frobenius operator by means of a 
Garlekin projection type approximation in order to compute the eigenvectors of the discretized 
operator corresponding to the eingenvalue 1. Although it can be used to approximate any SBR 
measure, the method is especially suited to approximate and study (almost) cyclical behavior 
of dynamical systems. However, its complexity makes the efficient implementation troublesome. 
A much simpler idea, which we shall adopt here, is to approximate the measure by truncating 



equation (2.1 ) for a reasonably large value of n. That is, we consider the approximating measure 

^ n— 1 

IJ.n{A;s,xo) ^ -^6T^k(^xo)(A) (5.1) 



n 

k=0 



which converges in a weak sense to as n tends to infinity, for almost all initial points xq G I 
and all /i^-continuity sets A. The iterations of Tg are known to be unstable with respect to the 
initial point in the sense that, given a small e > and a point x G (0, 1), the trajectories T^{x) 
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and Tg{x + e) become far apart exponentially fast. The approximation (5.1), however, is quite 
stable with respect to the initial point xq for large n. For instance, in Figure 5.1 we show the 
measure of the sets [0.1,0.2] and [0.4,0.6] obtained by using /i„(-;s,xo) with s = 0.5, for 50 
different initial points xq and 3 different truncation points n G {300, 000; 1, 000, 000; 3, 000, 000}. 
All plots are in the same scale (within set) in order to make comparison possible. In Table [5^1] we 
show basic statistics related to Figure [STTj Notice that, in average, the 1,000,000 and 3,000,000 
iteration cases are very similar and all cases are fairly stable with respect to the initial points 
(observe the scale). 



Set [0.1,0.2] with 300,000 Iterations 



Set [0.4,0.6] with 300,000 iterations 




02 04 



Set [0.1,0.2] with 1,000,000 iterations 




04 Oj6 



Set [0.1,0.2] with 3,000,000 iterations 



n 










QJQ Q2 0.4 06 






Set ^>l,0.e] wi#i 1 ,{H>D,{HH) lt«-aff<m& 




- '-'vo^i' j^v' Vv\ ' 


V 




OjQ Qi 0.1 Oj6 




ra 


Set [{].4,Q.fi.l withS.a'Qfl.a'Ofl Iterations 











oja CL2 a.4 Of a.s i.a 

tampon 



Figure 5.1: Performance of the approximation ( |5.1[ ) for truncation points n £ 
{300,000; 1,000, 000; 3, 000, 000} (top, middle and bottom, respectively) and 50 different initial 
points for s = 0.5. The measured sets are (left) [0.1,0.2] and (right) [0.4,0.6]. All plots within the same 
set are in the same scale. 



Table 5.1: Summary statistics for the data presented in Figure 5.1 



Set 


n 


300,000 


1,000,000 


3,000,000 


CO 


[min, max] 


[0.12511,0.13067] 


[0.12431,0.12901] 


[0.12688,0.12825] 




range 


0.00556 


0.00470 


0.00137 


,d. 


mean 


0.12790 


0.12775 


0.12777 


to" 


[min, max] 


[0.15349,0.16092] 


[0.15326,0.15944] 


[0.15676,0.15857] 


o 


range 


0.00743 


0.00618 


0.00181 


,d. 


mean 


0.15792 


0.15771 


0.15771 



Next question is how good is the approximation ( |5.l[ )? One way to test this is by test- 
ing whether the approximation is invariant under Tg. For given initial points, say xi, • • • ,Xk 
and some interval [a, 6], we calculate //„([a, 6]; s, Xj) and ^„(T~^([a, 6]); s, Xj). If the differ- 
ence between the two quantities is small for different pairs 
approximation is reasonably good. In Table 



1X7', Xn 



one can conclude that the 



fin{Ts ^{[a,b]);s,Xj) 



5.2 



point was taken to be 3,000,000 and s 



we present the difference 
for 7 different initial points and 3 different sets [a, b] 

0.5. From Table 5.2 we conclude that the approxi 



fin{[a,b];s,Xi) - 
The truncation 
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mation (5.1) performs very well in all cases and that it can be taken to be T^-invariant. As 
expected, when the differences are the smallest (< 10 ® in all cases). 

Table 5.2: Difference |/i„([a, 6]; s, Xi) — /i„(T^^^([a, &]); s, Xj) | for different values of xq and sets [a,b]. 
The truncation point was taken to be n = 3, 000, 000 and s — 0.5. The initial points are {xi, • • • , xy) — 
(tt, tt/{V2 + 1), TT A TT + A %/7, TT + 77, + ^/7) (mod 1). 





initial 


XI 


^2 


^3 


X4 


^5 


xe 


XJ 




















[0.05,0.2] 


XI 

3=3 

Xi 

3=5 

XQ 

X-j 


0.00000 
0.00019 
0.00040 
0.00008 
0.00004 
0.00062 
0.00022 


0.00019 
0.00000 
0.00000 
0.00030 
0.00020 
0.00043 
0.0004 


0.00040 
0.00020 
0.00000 
0.00047 
0.00044 
0.00022 
0.00062 


0.00008 
0.00027 
0.00047 
0.00000 
0.00003 
0.00070 
0.00015 


0.00004 
0.00024 
0.00044 
0.00003 
0.00000 
0.00066 
0.00018 


0.00062 
0.00043 
0.00022 
0.00070 
0.00066 
0.00000 
0.00084 


0.00022 
0.00042 
0.00062 
0.00015 
0.00018 
0.00084 
0.00000 






















initial 


XI 


X2 


^3 


XjX 


^5 


XQ 


x-r 




















[0.3,0.8] 


XX 
^2 
^3 
X4 

X5 

xe 

X-j 


0.00000 
0.00019 
0.00011 
0.00009 
0.00052 
0.00036 
0.00155 


0.00019 
0.00000 
0.00008 
0.00028 
0.00033 
0.00016 
0.00136 


0.00011 
0.00008 
0.00000 
0.00020 
0.00041 
0.00024 
0.00144 


0.00009 
0.00028 
0.00020 
0.00000 
0.00061 
0.00045 
0.00164 


0.00052 
0.00033 
0.00041 
0.00061 
0.00000 
0.00016 
0.00103 


0.00036 
0.00016 
0.00024 
0.00045 
0.00016 
0.00000 
0.00119 


0.00155 
0.00136 
0.00144 
0.00164 
0.00103 
0.00119 
0.00000 




initial 


XI 


^2 


^3 


X4 


^5 


xe 


XY 




















[0.7,0.95] 


XX 
3=2 
3=3 
x^ 

3=5 
XQ 
XJ 


0.00000 
0.00011 
0.00005 
0.00012 
0.00003 
0.00012 
0.00089 


0.00011 
0.00000 
0.00016 
0.00022 
0.00013 
0.00022 
0.00078 


0.00005 
0.00016 
0.00000 
0.00006 
0.00003 
0.00006 
0.00094 


0.00012 
0.00022 
0.00006 
0.00000 
0.00009 
0.00000 
0.00100 


0.00003 
0.00013 
0.00003 
0.00009 
0.00000 
0.00009 
0.00091 


0.00012 
0.00022 
0.00006 
0.00000 
0.00009 
0.00000 
0.00101 


0.00089 
0.00078 
0.00094 
0.00100 
0.00091 
0.00101 
0.00000 



In the remaining of this section we shall assume that s G (0, 1) has been fixed and xq £ (0, 1) 
has been chosen so that the approximation (5.1) converges to fis- Since no confusion will arise, 
we shall drop s and xq from the notation and write the approximation (5.1), based on a size n 
iteration vector, just by iJ^ni-)- 



Approximating ^ and the nodes of T, 



In order to approximate Fq ^, one can use an empirical version based on the same iteration 
vector from which /u„ is derived. First we need to define an approximation to Fq from which 
an approximation to Fq^ will be derived. Let Fn be the empirical distribution based on a 
size n iteration vector (xq, T's(xo), • • • ,T"~^(xo)) and let xi,-- - ,x„ be the jump point^ of 

Fn- Consider the set L„ := {0 = xo,xi,--- ,x„,x„+i = 1}. Given x G I\Ln, there exists a 
k £ {0, - ■ ■ , n} such that x G (xfc, Xfc+i). We define the approximate value of Fq{x), denoted by 
Fn{x), as the linear interpolation of x between the points (xfc,F„(xfc)) and (x^+i, F„(xa;+i)) , 
that is, we set 

p , ^ I Fn{xk+l) - Fnixk) \ ^ Fn{xk)Xk+l ~ Fn{xk+l)xk 

\ Xk+1 - Xk J Xk+1 - Xk ' 

If X G Ln, we simply define Fn{x) := Fn{x). Notice that, for each n, F„ :/—)•/ is a one-to-one, 
increasing and uniformly continuous function, so that its inverse, F'"^, is well defined and is also 
one-to-one and uniformly continuous. In the next proposition, we show that Fn{x) — t- Fq[x) and 
F~^{x) — )• Fq~^(x), both limits being uniform in x. 



^by the choice of xo, there wiU be exactly n jump points. 
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Proposition 5.1. Let Fn be the empirical distribution based on an iteration vector (xo,Ts(xo), 
• • • , r"~^(xo)) and let xi, - ■ ■ ,Xn be the jump points of Fn- Let Fn be the approximation ( |5.2[ ) 
based on {xi, ■ ■ ■ ,Xn} and F~^ be its inverse. Then, 

F„{x)^Fo{x) and F'^ix) ^ F^'ix), 

uniformly in x. 

Proof: By the Glivenko-Cantelli theorem, Fn{x) — )• Fq{x) uniformly in x G [0, 1], so that, given 
e > 0, one can find tiq := no(e) > such that if n > no, then — -Fo(2;)| < £ uniformly in 

X. Now, for x G (0, 1) (if x equals or 1, the result is trivial), there exists a A; G {1, • • • , n} such 
that x G [xfc, Xfc+i). Hence, if n > no 

|F„(x)-Fo(x)| < \F,,{x)~F,,{x)\ + \Fn{x)^Fo{x)\ 

< \Fnixk+l) ~ Fn{xk)\ + e 

< sup {\Fn{Xi+i) - Fn{Xi)\} + e 



< 



I,"' ,n— 1 

1 

n 



uniformly in x. To show the convergence of the inverse, let y G [0, 1] and e > be given and 
notice that F~^ being uniformly continuous, one can find a 6 := 5(e) > such that 

\x-y\<5 =^\F-\x)-F-\y)\<e. 

Now, since -F„ converges uniformly to Fq, there exists ni := ni{e) > such that, 

n> ni =^ \Fn{x) - Fq{x)\ < 5, 

for all X £ I. Also, since Fq is one to one, there exists vq G [0, 1] such that y = Fq{vq). Therefore, 
if n > ni 

\F-\y) - F^\y)\ = \F-\Fo{vo)) ~ vo\ = \F-\Fo{vo)) - F-\FM)\ < e 

and since ni is independent of y, the desired convergence follows. ■ 

As for the end points {ah,k}'kLo of the nodes of T^, let {xi,-- - ,Xm} £ (0,1), Xi ^ xj 
and consider the set {T^{xi), ■ ■ ■ ,T^{xm)}, for m > sufficiently larg^ Note that Ohfl = 
and af^2h = 1, for any h. Let D = {i : T^{xi) > T^^(xj+i)} C {!,••• ,rn}. The set D 
contains the indexes i G {1, • • • , m} for which the interval [xj, Xj+i] contains a discontinuity of 
Tg . Let denote the ordered elements of D, so that the interval [x^^. ,Xd^.^J contains 

the j-th discontinuity of T^. Now consider the function T*^.^ : [x^. , x^^^+i] — )• [0,2] given by 
T*j^{x; s) := Tg~^{x) + (Tg"^ {x))^~^^ and notice that we can write 

T^{x) = T*,^{x;s)-S[^^,]{T*^ix:,s)). 

Since there is a discontinuity of in the interval we have T*^(xdjs) < 1 and 

'^ihi^di+i', s) > 1 and since T*^ is continuous and increasing, there exists a point x G [xd^,Xdi+i] 



^By "sufficiently large" we mean that m should be at least large enough to guarantee that the set 
{T^{xi),--- ,T!^{x,n)} reflects the 2'' — 1 discontinuities of , or, in other words, m > 2''. The limits in m 
taken for an approximation are understood to be in terms of partitions, that is, we start with a sufficiently large 
set of points, say Im = {xi,--- ,Xm} and consider refinements of the form Im+i ~ Im[J {xm+i}, ■ ■ ■ ,Im+k = 
Im+k-i[J {xm+k}- Suppose that R,n := R{Im) is an approximation based on Im- For a sequence of refinements 
I ]^ we consider the sequence {ii(/fe)}^„, J. Whenever the last limit exists, we set lim Rm ~ lim R(Ik)- 

m — ^oc k — yoo 
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such that T*{x;s) = 1, which is precisely ah^i- With this in mind, let a™^ denote the approx- 
imation to ah,i obtained from {xi, • • • , Xm} by using a linear interpolation between the points 
{^d,,T*hixd,;s)) and {xd^+i,T*f^{xd,+i; s)) . That is, a^- is given by 

= + , """'^I'l':. ^ (1 - T*,ix,^;s)), (5.3) 

for all di £ D. Clearly aV^- — > Uhi, since — XdA — > and by the continuity of T*,, 

for each iE {!,••• ,2^-1}. 



Approximating Th,k 

Concerning the approximation of Th,k, we shall use an argument based on an empirical inverse 
and linear interpolation, but we shall also need a doubling argument in order to improve accuracy 
of the approximation near the discontinuities and guarantee the uniform convergence of the 
approximation to its target. So let {0 = xi, • • • ,Xm = 1} £ /, < Xj and consider the set 
[Tg{xi), ■ ■ ■ ,T!!'{xm)}, for m > sufficiently large. Given y G [0, 1], recall that the inverse of 
Tg{y) is a size 2^ vector which we denoted by (Thfl, ■ ■ ■ >7/j 2'^-i)- Let again D = : Tg{xi) > 
Tg{xi^i)^ C {1, • • • , m} and {dj}^^^^ be the ordered points in D. Suppose that we know exactly 
or have good estimates for the nodes {ah,k}tLo of (fo'^ instance, we could use {a^^}|!lQ, as 
described before, based on the same set {xi, • • • , Xm} considered here). For i = 0, • • • , 2^* — 1, 
let 

and 

Given y € [0, 1], for each i = 0, • • • , 2^^ — 1, there exists a y\^l e I^^i such that y € [yl^lTvl^^^^]- 
We define the approximation Tfi^^{y) of Th,i{y), as being the linear interpolation of y between 

the points [x^^J , y^j^J) and {x\^J'^j Vh/^^) ■ That is, for each i = 0, • • • , 2^" - 1, 

(fe+i) _ (fe) 

Vh.i Vh.i 

Notice that if y equals or 1, we have Tf^^{y) = Th,i{y)- Also, as the partition {xi, • • • ,Xm} 
increases, Ix^+i — x^l — > and the uniform continuity of T!^ clearly implies T^Ay) — > 

Th,i{y), for each y G [0, 1], for i = 0, • • • , 2^* — 1. More is true: the convergence is actually 
uniform in y, as we show in the next proposition. 



Proposition 5.2. Let 7^™ be the approximation ofTh^k given by (5.4) based on a partition Rm- 
Then, 

foreachk = 0, • • • ,2^^-!, as m goes to infinity (that is, as the partition gets thinner). Moreover, 
the convergence is uniform in y G [0, 1]. 

Proof: Given e > 0, the uniform continuity of 7h,k implies the existence of a 5 := 5(e) > 
such that 

\x-y\ < S =^ \Th,k{x) ~ Th,k{y)\ < e, 
for all X G [0, 1]. Let i?o = {0 = xi, • • • , Xmo = 1} G / for a sufficiently large mo G M* such that 

sup {\x.,+i - Xi\} < S. 

- .mn — 1 
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For m > mo, let Rm = {x*, • • • , x^} D Rq he a size m refinement of Rq. Given y £ (0, 1), for 
each i = 0, • • • , 2^ — 1, let 7^"J be the approximation (5.4) based on Rm- By construction and 
since y G (0, 1), it follows that 

%4xi'}) < T^^{y) < r..«+^)) and < %m < TuMr^)^ 

so that 

\Tay)-rhAy)\ < IThMT^) ~ r.M'.Dl 

< sup {\Th,i{^j+i) -ThA{xj)\} < e, 

J — I,--- ,771—1 

for all y G (0,1). If y G {0)1}) by construction Th,i{y) = T^^{y), so that the result follows 
uniformly for all y G [0, 1], as desired. ■ 

5.1 Approximating the lag h MP copula 



With these approximations in hand, we can now define the approximation for the copula Cxt,x, 



when if is almost everywhere increasing given in Proposition 3.2 but in the form (3.5). For 
(tt, u) G n > and m > 2^, we set 



t+h 



^ fc=0 ^ 

+ A*77 ( , min {F-\u) , TCno {Fn\v)) }] ) , (5.5) 

where ng := no(m, n) = {fc : u G [-^ri(a™fc)' -^"('^/Tfc+i)) } ^™ ^0 ~ '^o since F„ converges 



uniformly to Fq and converges to a/i,fc- In the next theorem we establish the convergence of 



the approximation (5.5) to the true copula. 



Theorem 5.1. Let Cm,n{u,v; h) he given by (5.5). Then, for all {u,v) G , t>0 and h> 



lim lim Cm,n{uTv;h)= lim lim Cm,n{u,v'ih) = lim Cm,n{uT'V]h) 

71— VOO 777— >-00 ' 777— >-00 77— >00 ' 777 , 77 —f CO ' 



and the common limit is Cxt,Xt+h{u.,v) (given by (3.4)j. Furthermore, the limits above are 
uniform in (n, f) G 



The proof of Theorem 5.1 , is a consequence of the following stronger lemma. 



Lemma 5.1. Let {/XnjrieiN be a sequence of probability measures defined in I such that 
/i. Let fn : I ^ I be a sequence of continuous functions converging uniformly to a function 
f : L ^ L. Let {am}mefi be a sequence of real numbers such that Um G [0, 1] for all m and 
ttm — ^ o- Also let gm ■ [am, 1] ^ L be a sequence of continuous functions converging uniformly 
to a function g : L ^ L , Sm,n{v) := [um, gm{fn{v))] and S{v) := [a,g{f{v))\. Then, 

lim lim /i„(5™.„(u)) = lim lim ^„(S'„,„(i;)) = lim /i„(5™.„(u)) = ^(S'(w)) 

777— >00 77— fOO ' 77— fOO 777— >-00 777,77—^00 \ ■ / 

uniformly in v £ L. 

Proof: For all m,n > and v G [0, 1], let Sm,n{v) and S{v) be as in the enunciate and let 

Sn{v) := [a,g[fn{v))] and ^^(t;) := [a„,5„(/(u))] • 

Notice that all sets just defined are //-continuity sets for all m, n and v. Since the convergence 
of /„ to / is uniform, we have 

lim g,Jf^{v))^ lim lim gm(fn{v))= lim lim g„Jfn{v)) = g(f{v)) 



m.n~^oo n-^oo m^oo m— ^oo n— ^oo 
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for all V, so that, both, the iterated and the double limits exist and Sm,n{v) — )• S{v), for all 
V G [0, 1]. Also notice that we have Ss^ni^) — uniformly in m, n and and since 

converges weakly to ^ and / is a //-continuity set, it follows that 

j 5i{x)diJLn — > j 5i{x)dii. 

Now, in one hand, since Sm,n{u) Sjn{v) for all v and Ssm,n ^ ^i^^Y Lebesgue convergence 
theorem, it follows that 

and, since < and J (J/d/x < oo, by the Lebesgue dominated theorem, we conclude that 

j ^^rnix)d^. -^^ j Ss{x)dfi = iJ,{S{v)), 



which shows that lim lim l-Ln(Smn{v)) = fi(S(v)) and the convergence holds uniformly in 

V G (0, 1). On the other hand, since ds^„ < and J 5/d/x„ < oo, by the Lebesgue dominated 
theorem, it follows that 



(v)) = / 6s^^{x)diin — > / Ss„{x)dnn, 

J ' m-)-cx> J 



and, since < 5/ and / Sjd^n J ^idfi, by the Lebesgue convergence theorem we conclude 
that. 



j 5s„ (a;)d/i„ j Ss{x)diJ, = iJ,{S{v)), 



that is, lim lim fin(Smn{v)) = l^iS{v)), which also holds uniformly in v. Since the iterated 

limits are established, in order to finish the proof we need to show that the double limit exists 

and is equal to the iterated ones. Let e > be given. Since fi <^ \, the Radon-Nikodym theorem 
implies the existence of a non-negative continuous function h, which will be bounded since we 
are restricted to the interval /, such that, for any A G B{I), 



At(A) = j h{x)dX<MX{A), 



where M = sup{h{x)} < oo. Now, since am — >■ a, one can find mi := mi{e) > such that, if 
m > mi, 

r s e 



and 



M(Jri(e)) <ma( 



lOM' 



lOM 

£ 



The uniform convergence of gm to g implies the existence of m2 '■= m2(e) > such that, if 
m > m-2, \gm{x) — g{x)\ < £/20M, for all a; G /, or equivalently, taking x = fn{v), if m > m2 



(/n (■'')) £ 5 (/n (■'')) 



..9(/n(tO) 



^ 20MJ 

Now, the uniform continuity of g implies the existence ol a 6 := 6{e) > such that 

£ 



\x- fn{v)\ < 5 



\9{x) - g{fn{v)) \ < 



20M 



But since fn converges to / uniformly, there exists a ni — ni{5) > such that 



n > ni 



\fn{v) - f{v)\ < 6, 
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for all V so that, taking x = f{v), for n > ni, we have 



9{fn{v))e[g{f(v))-—,g{f{v)) 
for all V G I. Hence, if we take m > m2 and n > ni, 



20M 



Aln{v)) 



20M 



and 

so that, setting 



g{f{v)),g{f{v)) 



e 

lOM 



-?f2(£) \g(.f(v)) —,g(fiv)) + 



for m > m2 and n > ni, it follows that 



for all V G I. Also observe that 



KMe)) < Mx{[g{fiv)) - ^,g{fiv)) + 



lOM 



< 



The convergence of fin to /u implies the existence of n2 := ^^2(£) > such that if n > n2 (i^j(e) 
is a /i— continuity set) 

\fj.„{K,{e)) - fi{K,{e)) \ < |, 

for i = 1, 2. Also, if we set Fn{x) = /u„([0, x]) and -Fo(2;) = /u([0, x]) , then Fq is continuous (since 
/U ^ A), Fn — ?• -Fo, and, by Polya's theorem, there exists a ns := n3(e) > such that, if n > na 

sup{|F„(a;)-Fo(a;)|} < ^ 



10 



Now, notice that, if n > 723 

\fi4S{v)) - fi{S{v))\ < \F4g{f{v))) ^ F„{g{f{v)))\ + \F„{a) ~ Fo{a)\ 
< 2sup{|F„(x)-Fo(x)|} < ^, 



xei 



for all V G L Observe further that, by construction, if m > max{mi,m2} and n > rii, 

SmAv)\Siv) d Ki{e)ljK2{e), 
for all V so that, setting no = no{e) := max{mi, m2, ni, 71,2, na}, if m,n > no, we have 

|Mn(5'm, -^(S'(w))| < 

< \fln{S„i,n{v)) - fln{S{v)) \ + |/Z„(S'(w)) - fl{S{v))\ 

< \fln{Kl(e)) + fln{K2{e))\ + I 



< 
< s 



fin{Ki{e)) -fi{Ki{e)) \ +fi{K^{e)) (^2(e)) + |M(^2(e)) (^2(e)) | + ! 



for all 11, which implies the existence of the double limit, equality with the iterated ones and the 
desired uniform convergence. ■ 



Proof of Theorem 5.1: First notice that taking fn = F^^,gm= Tf^f., am = a^^^, it follows 
from Lemma Is. II that 

l^n[[al,,:r^k{F-\v))]) ^-^^ /i([a,,fe,7;a-(i^o"' («))]), 
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for each k = 0, - ■ ■ , tiq — 1. It remains to show that 

^lim^ fin ( [aZn^ , min {F^i rr„o {F-\v)) }] ) - , min {F„-' {u),Th. no {Fo'\v)) }] 



and that the iterated limits exist and are equal to the double limit. First, since we can write 
min{ti,t;} = — ^"^"^ ; it is routine to show that if /« — >■ / uniformly, with and / 
uniformly continuous and gm 9 uniformly, with Qm and g uniformly continuous, we have 
min I /„(«), (/n(^^)) } converging uniformly to min {/(u), 17 (/(f)) | in n, m, u and v. So, 
the problem simplifies to show that if am — >• a, gm,n{u,v) is a sequence of functions such that 

gm,n{u,v) g{u,v) Uniformly in u,v,n,m and am < gm,n{u,v) for all u,v,n,m and /x„ //, 
then 

lim fin([am,gm.n{u,v)]) ^ fi([a, g{u, v)]) , 

m,n— fcjo 

uniformly in u and w and the double limit above is equal to the iterated limits. A similar 
argument to the one used in Lemma |5.1| to establish the existence and equality of the iterated 
limits can be used to show the existence and equality of the iterated limits in this case. As 



for the double limit, let M be as in the proof of Lemma 5.1 By the uniform convergence of 
9m,n{u^v) to g{u,v) and since gm,n and g are uniformly continuous for all m,n, it follows that 
there exists mi := mi{e) > 0, depending on e only, such that, if m,n > mi, 

.gm.nlw,-") e K{e) := g{u,v) ~ j^,giu,v) + 



for all u and v and fi(^K(e)^ < e/5. The rest of the proof is carried out by mimicking the 



proof of Lemma 5.1 with the obvious adaptations. Identification of gm,niu, v), g{u, v), am and a 
with min {F-i(n),7;7„^(F-i (-(;))}, mm{FQ^{u),Th,no{Fo^{v))}, a^^^ and ah,no, respectively, 
completes the proof. ■ 



Remark 5.1. Notice that neither the convergence proved in Lemma 5.1 nor the one in Theorem 
15.11 is uniform in m and n. 



As for the case when ip is almost everywhere decreasing, we observe that, in view of (3.7), 
the function 

Cm.„(", v; h) = u + V - I + Cm,„(l -u,l- v;h) 



is an approximation to the copula in (3.6). Clearly Cj^„ converges to the true copula as m and 
n tends to infinity (view either as an iterated or a double limit) and the convergence is uniform 
in {u, v). 



Implementation and Random Variate Generation 

The implementation of the approximations so far discussed is routine. All the approximations 
we mentioned can share the same iteration vector, which further improves the efficiency and 



precision of the task and greatly reduces the computational burden. In the top panel of Figure 5.2 



we show the three dimensional plot of the lag 1 and 2 MP copula for values of s G {0.1, 0.4}. The 



respective level plots are shown in the bottom panel of Figure 5.2 Notice the non-exchangeability 
of the copulas in all cases. 



Obtaining random samples from an MP copulas is a trivial task in view of Proposition 3.4 
There we show that the support of an MP copula is the union of graphs of certain linear functions. 
The following algorithm can be used to generate a pair of variates from a bidimensional MP 
copula for ip an almost everywhere increasing function. 



1. Generate an uniform (0, 1) variate u. 
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Figure 5.2: From left to right, three dimensional plots of the lag 1 MP copula for s £ {0.1,0.4} and lag 2 MP 
copula for the same parameters (top panel) and respective level sets (bottom panel) obtained from approximation 



2. Let kq denote the index for which u G [Fo(a/i_Ko)) -^o(a/i,Ko+i)] and set v = ^^^^(w)- 

3. The desired pair is {u,v). 



In practice the T^-invariant probabiHty measure is unknown and Fq has to be approximated. 
Furthermore, most of times the nodes related to T^, for /i > 0, s E (0, 1) cannot be analytically 
obtained. However, we can apply the approximations developed in this section together with 
the algorithm above to obtain approximated samples from MP copulas. In Figure [53| we show 
500 approximated sample points from a lag 1 and 2 MP copula for s G {0.1,0.4} and if an 
almost everywhere. Obvious modifications in the algorithm, allow handling the case where (p is 
an almost everywhere decreasing function. 



500 sample points of the lag 1 MP copula for s 



500 sample points of the lag 1 MP copula for s 





0.0 0.2 0.4 0.6 0.8 1.0 



0.0 0.2 0.4 O.S 0.8 1.0 



500 sample points of the lag 2 MP copula for s = 

:] / / / / 

/ / / / 

- / / / / 
'Is / 

= ' ! 

= - f s 

0.0 0.2 0.4 0.6 0.8 1.0 



500 sample points of the lag 2 MP copula for s = 0.4 



0.0 0.2 0.4 O.S 0.8 1.0 



Figure 5.3: Left to right: 500 approximated sample points from a lag 1 MP copula for s € {0.1,0.4} and lag 2 
MP copula for the same parameters. 



Remark 5.2. For small values of the lag, the resemblance of the sample to a piecewise contin- 
uous function is very clear, but this is not always the case as it can be seen in Figure 5.4, where 
we show 500 approximated sample points of the lag 4, 5 and 7 MP copulas for s = 0.2. This 
is a general principle, for a fixed sample size the higher the lag, the harder to distinguish the 
support of the copula based on the sample, since the number of branches of grow as fast as 
2^. For instance, for /i = 7 in Figure 5.4 is difficult to say that the sample came from a singular 
copula at all. 



S.R.C. Lopes and G. Pumi 



21 



500 sample points of the lag4MPcopijlafors = 0.2 



500 sample points of the iag 5 MP copuia for s = 0.2 



500 sample points of tiie iag 7 MP copuia for s = 0.2 



0.0 0.2 0.4 0.6 0.8 1.0 



0.0 0.2 0.4 0.6 0.8 1.0 



0.0 0.2 0.4 0.6 0.8 1.0 



Figure 5.4: Left to right: 500 approximated sample points from the lag 4, 5 and 7 MP copulas for s — 0.2. 



6 Application 



In this section we apply the theory developed in Section 3 to the problem of estimating the 
parameter s in MP processes. This problem have been studied before in Olbermann et. al 
(2007), where the authors adapt and apply several estimation methods from the classical theory 
of long-range dependence to the problem of estimating the parameter s. In this section we 
propose an estimator for the parameter s based on the ideas developed in Section 3, which is 
both, precise and fast. 

The mathematical framework is as follows. Let s £ (0, 1) and consider the associated MP 
process {Xn}neJN for 93 the identity map. Suppose we observe a realization xi, - ■ ■ ,xn from Xn 



and our goal is to estimate the unknown parameter s. Let a := a{s) G 



'1 v^-l 

.2' 



discontinuity point of the MP transformation and notice that s and a are related by 



) denote the 



a + a 



l+s _ 



log(l - a) 
log(a) 



1. 



Hence, the problem of estimating s is equivalent to the problem of estimating a. 



To define the proposed estimator, we start by observing that Proposition 3.4 for h = 1 
implies that the lag 1 MP copula's support is given by the graph of the piecewise linear function 



Fo(a) ' 
x—Fo{a) 
l--Fo(a) 



if [0,Fo(a)) 

if [Fo{a),l]. 



so that, any (independent or correlated) sample from a lag 1 MP copula consists of points 
scattered through the lines defined by i (see Figure 5.3 ). The discontinuity point of the function 
is precisely Fq (a). Let 



Fo{xi), for z = 1, • • • , A^, and consider the series {ui := (ui, Ui+i)} 



N-l 
1=1 ■ 



By Sklar's Theorem, {uj}^^^ is a (correlated) sample from the lag 1 MP copula, so all points 
should lie in the graph of the function i. 

These considerations suggest the following procedure to obtain s based on a path xi, • • • , xat 
of Xn within a given accuracy e > 0. We choose sq G (0,1) as an initial guess for s and calculate 
{xf, So), i = 1, • • • ,N, where Fn is the approximation of Fq given in ( 5.2 ). Next we define 
from which we estimate the slope of the two branches of the approximated 



sample from the lag 1 MP copula obtained by this way. The discontinuity point (and hence s) 
can then be easily calculated. In this manner we obtain an estimative s which can be compared 
to Sq. If So is close to the true value s, then the difference between s and sq should be small. If 
not, we choose another starting value and repeat the operation until obtain the desired accuracy. 
This leads to an optimization procedure to obtain s within a predefined accuracy. 



To illustrate the procedure. Figure |6.1[ a) shows a sample path of an MP process for s = 0.2, 
with N = 200 while Figure |6.1[ b) shows the sample path yi = Fn{xi; 0.2), i = 1, • • • ,N. From 
{yi}iLi, we construct the sequence where m = [Fn{yi] s) , Fn{yi+i; s)) , i = !,■■ ■ ,N- 
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1, for the correctly specified s = 0.2 and for s = 0.3. Figure pTTt c) presents the graph of 
{uj}^^^ obtained from the correct specification of s, while Figure p. iK d) shows the graph of the 
misspecified one. In Figures 6.1 ^c) and 6.1 d), the solid lines represent the respective theoretical 

Some distortion in the points can be seen given to 



support of the copula given in Proposition [3^ 
the use of t he approximation Fn instead of the theoretical Fq, especially in lower quantiles. From 
Figure 6.1 'd) it is clear that the line obtained from the sequence {ui}^is[^ and the theoretical 
one for the chosen value of sq, namely, 0.3, do not match, while for the correct specified one in 
Figure |6?T|c), they do. 




Figure 6.1: (a) Sample path 2:1, •• • , 2:200 of an MP process with s — 0.2 starting at v^(modl). (b) Sample path 
Hi ~ Fn{xi). Plot of Ui = {yi,yi+i) for the (c) correct and (d) misspecified s. The solid lines correspond to the 
theoretical support of the respective lag 1 MP copula. 



The procedure just outlined is, however, computationally expensive given the fact that to 
calculate the approximation Fn with reasonable stability and accuracy, for each s, it requires 



the construction of an iteration vector of large size (see Figure 5.1 and Table 5.1 ). Such an opti- 
mization procedure can easily take hundreds of evaluations, depending on the desired accuracy, 
and hence, can be a very time consuming task. 



To overcome this difficulty, observe that in Figures |6.1[ a) and |6.1[ b) , little differences can 
be seen between them. In fact, since Fq is a smooth distribution, an alternative is to apply 
the previous argument to the points Vi := {xi,Xi+i), i = 1, • • • ,N. There will certainly be 
some distortion in the lines due to the absence of Fq, but we expect to be able to estimate the 
discontinuity point a based on Vi by similar idea as before. 



As an illustration. Figure 6.2 shows the plots of Vi = (xj, Xj+i), i = 1, - ■ ■ , 199, based on MP 
processes with s € {0.2, 0.4, 0.6, 0.8} all starting at \/5(modl). The solid lines correspond to the 
lines joining the points (0, 0) and (a, 1) and joining (a, 0) and (1, 1), where a denotes the correct 
discontinuity point of the respective MP transformation. From the graphs in Figure 6.2 we see 
the identification of the line based on Vi with the correct line, especially in the second branch 
of the graph. That is so because a G (j, ^~^ ) , so that the second branch, being smaller, is less 
affected by the distortion due to the absence of Fq. 

In order to assess the performance of the estimation procedure, we perform the following 
experiment. We randomly select 100 initial point^in (0, 1) and for each initial point we generate 
a path (of size N = 200) of an MP process for s £ {0.1, 0.15, •• • ,0.95}. For each path, say 
xi,- " 5 X200, we perform the proposed estimation procedure. In order to estimate a, we applied 
two methods: the first one is a simple least squares method applied to the points lying in 
the second branch of (xi,Xi+i). The second method is the following: let {xmo, Xmo+i) and 
{xmi,Xmi+i) denote the points among the ones lying on the second branch of {(a;^, Xj+i)}^"^"^ 



''Tables with the initial values applied in our experiments and the complete simulation results are available 
upon request. 
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Figure 6.2: Plot of Vi = {xi,Xi+i), i = 1, - ■ ■ , 199 from a sample path of an MP process with (a) s = 0.2, (b) 
s — 0.4, (c) s — 0.6 and (d) s — 0.8. The solid lines correspond to the lines joining the points (0, 0) and (a, 1) and 
joining (a, 0) and (1, 1), where a denotes the correct discontinuity point of the respective MP transformation. 



for which Xmo is minimum and Xmi is maximum. We define the estimator of a, say a, as 

a — — -, where A := rno+i B :— Xm^+i ~ Axmo- (6-6) 

For reference, in the subsequent we shall call this the min-max procedure. Geometrically, a is 
the inverse image of by the linear function joining {xmo,Xmo+i) and {xmi, Xmi+i)- 



Histogram of the L5 proc. for s = 0.5 



Histogram of the min-max proc. for s = 0.5 



(a) 



(b) 



(c) 



(d) 



Figure 6.3: Plot of the estimated values for s £ {0.1,0.5,0.9} for 100 random initial points by using (a) the 
least squares procedure and (b) the min-max procedure. The dashed lines correspond to the correct value of s. 
Also shown the histogram of the estimated values for s = 0.5 by using (c) the least squares procedure and (d) the 
min-max procedure. 



Table |6.1| summarizes the experiment results by presenting the mean, range, standard de- 
viation (st.d.) and mean square error (mse) of the results. Figures |6.3[a) and 



graphically the results for both methods for s S {0.1, 0.5, 0.9} while in Figures 6.3 
the histogram of the results for s 



6.3||b) present 
and |6.3[ d), 

0.5 are presented. From Table [6^1] and Figure [6131 we see 



that the min-max procedure (MM) outperforms the least squares estimates (LS) obtained. Some 
bias can be seen for both estimates, especially when s increases. 

The min-max procedure can be carried out even for time series of sample size as small as 20, 
as long as the second branch of {{xi, Xi-\-i)}^^S[^ contains at least 2 points, which does not always 
happen (for instance, for N = 110, a sample path of an MP processes with s = 0.8 starting at 
\/74(modl) has only one point in the second branch). In such a situation, a straightforward 
adaptation of the min-max procedure can be applied to the first branch and still yields reasonable 
estimates. The closer to and 1 the points Xmo and x^i in (6.6) are, respectively, the better 
the estimation performance. 



24 



Copulas Related to Manneville-Pomeau Processes 



Table 6.1: Summary statistics of the experiment results. Presented are the mean estimate (s), the range, 
standard deviation (st.d.) and the mean square error (mse) of the estimates. The min-max procedure is denoted 
by MM while LS denotes the least squares. 



Proc. 


s 


s 


range 


st.d. 


mse 


s 


s 


range 


st.d. 


mse 


i\ /ri\ /T 
MM 

LS 


0.10 


0.1008 
0.1087 


[0.1000, 0.1024] 
[0.1056,0.1128] 


O.OOOd 
0.0017 


0* 
0.0001 


0.55 


0.5581 
0.6036 


[0.5500, 0.5888] 
[0.5501,0.6287] 


0.0069 
0.0142 


0.0001 
0.0031 


MM 

LS 


0.15 


U.loiD 

0.1632 


[O.ioUU, u.ioy 
[0.1573,0.1710] 


O.OOio 
0.0026 


U 

0.0002 


0.60 


0.6091 

0.6545 


[O.OOOO, U.64oi] 

[0.6012, 0.7023] 


U.UU90 
0.0186 


0.0002 

0.0033 


MM 
LS 


0.20 


0.2023 
0.2179 


[0.2001, 0.2101] 
[0.2098,0.2315] 


0.0021 
0.0038 


0* 
0.0003 


0.65 


O.ddOO 
0.7089 


[0.6501, 0.6927] 
[0.6579,0.7584] 


0.0087 
0.0197 


0.0002 
0.0039 


MM 

LS 


0.25 


0.2534 
0.2730 


[0.2501, 0.2632] 
[0.2636, 0.2875] 


0.0027 
0.0049 


0* 
0.0006 


0.70 


0.7125 
0.7646 


[0.7001, 0.7726] 
[0.7038,0.8170] 


0.0119 
0.0226 


r\ r\r\r\o 

0.0003 
0.0047 


MM 

LS 


0.30 


U.oUoO 
0.3272 


[U.oUUU, U.olzoJ 
[0.3128,0.3410] 


U.UUzo 
0.0052 


U 

0.0008 


0.75 


U. 

0.8177 


[u. ( oUz, U.oUiyj 
[0.7505,0.8612] 


n m 1 n 
U.UilU 

0.0246 


0.0052 


MM 
LS 


0.35 


0.3544 
0.3835 


[0.3500, 0.3669] 
[0.3550, 0.4024] 


0.0039 
0.0078 


0* 
0.0012 


0.80 


0.8165 
0.8781 


[0.8001,0.8659] 
[0.8005, 0.9449] 


0.0131 
0.0277 


0.0004 
0.0069 


MM 
LS 


0.40 


0.4050 
0.4367 


[0.4000, 0.4214] 
[0.4082, 0.4570] 


0.0049 
0.0078 


0* 
0.0014 


0.85 


0.8677 
0.9307 


[0.8500,0.9428] 
[0.8507,1.0111] 


0.0151 
0.0297 


0.0005 
0.0074 


MM 
LS 


0.45 


0.4556 
0.4909 


[0.4501, 0.4703] 
[0.4702, 0.5174] 


0.0046 
0.0102 


0.0001 
0.0018 


0.90 


0.9172 
0.9774 


[0.9002, 0.9704] 
[0.9011, 1.0477] 


0.0141 
0.0306 


0.0005 
0.0069 


MM 
LS 


0.50 


0.5065 
0.5475 


[0.5001,0.5189] 
[0.5059, 0.5746] 


0.0051 
0.0111 


0.0001 
0.0024 


0.95 


0.9706 
1.0371 


[0.9500, 1.0517] 
[0.9500,1.1641] 


0.0189 
0.0384 


0.0008 
0.0090 



Note: 0* means that the mse is smaller than 5 x 10 ^ . 



7 Conclusions 

In this work we derive the copulas related to Manneville-Pomeau processes for almost everywhere 
monotonic functions ip. In the bidimcnsional case, wc find that the copulas of any random pair 
{XtiX^jf-h) depend only on the lag h and are singular. The support of the copulas is derived as 
well. 

As for the multidimensional case, when ip is increasing almost everywhere, the functional 
form of the copulas are very similar to the ones in derived in the bidimcnsional case. We conclude 
that the copulas of vectors {Xt^ ,XtJ and (/Jq. T"*""*' (Uo),- ■ ■ , r*"-*i (C/q)) are the same. When 
(p is decreasing almost everywhere, we find that the copulas of an n-dimensional random vector 
from an MP process can be deduced from the ones derived for the increasing case. 

The copulas derived here depend on the T^-invariant measure //^ which has no explicit 
formula. For the bidimcnsional case, we propose an approximation to the copula which is shown 
to converge uniformly to the true copula. Prom this approximation, we are able to present plots 
of the copulas for different parameters and lags and to present a simple algorithm to generate 
approximated samples from the copulas. Some simple numerical calculation are presented to 
test the steps of the approximation. To illustrate the usefulness of the theory, we derive a fast 
estimation procedure of the underlying parameter s in Manneville-Pomeau processes. 

Acknowledgements 

Silvia R.C. Lopes research was partially supported by CNPq-Brazil, by CAPES-Brazil, by INCT em 
Matemdtica and also by Pronex Probabilidade e Processos Estocdsticos - E-26/170. 008/2008 -APQl. 



S.R.C. Lopes and G. Pumi 



25 



Guilherme Pumi was partially supported by CAPES/Fulbright Grant BEX 2910/06-3 and by CNPq- 
Brazil. 

References 

Billingsley, P. (1999). Convergence of Probability Measures. 2nd ed., New York: John Wiley. MR1700749 

Chazottes, J.-R., P. Collet and B. Sclimitt (2005). "Statistical Consequences of the Devroye Inequality 
for Processes. Applications to a Class of Non-Uniformly Hyperbolic Dynamical Systems" . Nonlin- 
earity, 18(5), 2341-2364. MR2165706 

Dellnitz, M. and Junge, O. (1999). "On the Approximation of Complicated Dynamical Behavior". 
SIAM Journal of Numerical Analysis, Vol. 36(2), 491-515. MR1668207 

Fisher, A.M. and Lopes, A. (2001). "Exact Bounds for the Polynomial Decay of Correlation, 1/f 
Noise and the CLT for the Equilibrium State of a Non-Holder Potential" . Nonlinearity, Vol. 14, 
1071-1104. MR1862813 

Joe, H. (1997). Multivariate Models and Dependence Concepts. Monographs on Statistics and Applied 
Probability, 73. London: Chapman & Hall. MR1462613 

Kolesarova, A.; Mesiar, R.; Sempi, C. (2008). "Measure-preserving transformations, copulae and com- 
patibility". Mediterranean Journal of Mathematics, 5(3), 325-339. MR2465579 

Lopes, A. and Lopes, S.R.C. (1998). "Parametric Estimation and Spectral Analysis of Piecewise Linear 
Maps of the Interval". Advances in Applied Probability, Vol. 30, 757-776. MR1663557 

Maes, C; Redig, F.; Takens, F.; Moffaert, A. and Verbitski, E. (2000). "Intermittency and Weak Gibbs 
States". JVonJinearitj, Vol. 13, 1681-1698. MR1781814 

Nelsen, R.B. (2006). An Introduction to Copulas. New York: Springer- Verlag. 2nd Edition. MR2197664 

Olberman, B.P; Lopes, S.R.C and Lopes, A.O. (2007). "Parameter Estimation in Manneville-Pomeau 
Processes" . Unpublished manuscript. Source: |arXiv:0707.1600[ 

Pianigiani, G. (1980). "First Return Map and Invariant Measures". Israel Journal of Mathematics, Vol. 
35, 32-48. MR0576460 

PoUicott, M. and Yuri, M. (1998). Dynamical Systems and Ergodic Theory. Cambridge: Cambridge 
University Press. MR1627681 

Royden, H.L. (1988). Real Analysis. New York: Macmillan. 3nd Edition. MR1013117 

Schweizer, B. and Sklar, A. (2005). Probabilistic Metric Spaces. Mincola: Dover Publications. MR0790314 

Young, L-S. (1999). "Recurrence Times and Rates of Mixing". Israel Journal of Mathematics, Vol. 
110, 153-188. MR1750438 

Zebrowsky, J.J. (2001). "Intermittency in Human Heart Rate Variability". Acta Physica Polonica B, 
Vol. 32, 1531-1540. 



